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ABSTRACT 

The Farley-Buneman instability is studied in the partially ionized plasma 
of the solar chromosphere taking into account the finite magnetization of the 
ions and Coulomb collisions. We obtain the threshold value for the relative 
velocity between ions and electrons necessary for the instability to develop. It 
is shown that Coulomb collisions play a destabilizing role in the sense that they 
enable the instability even in the regions where the ion magnetization is greater 
than unity. By applying these results to chromospheric conditions, we show 
that the Farley-Buneman instability can not be responsible for the quasi-steady 
heating of the solar chromosphere. However, in the presence of strong cross-field 
currents it can produce small-scale, ~ 0.1 — 3 m, density irregularities in the solar 
chromosphere. These irregularities can cause scintillations of radio waves with 
similar wave lengths and provide a tool for remote chromospheric sensing. 

Subject headings: Sun: atmospheric motions — Sun: chromosphere 
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Introduction 



The mechanism of chromospheric heating is a major puzzle of solar physics since it was 
discovered that the temperature in the solar chromosphere is much higher than what can 



be expected for a plasma in radiative equili 
of chromospheric heating was proposed by 



Drium. The 



irst possib 



Biermann 



(119461 ) and 



e scenario for explain ing 



Schwarzschild 



(119481 ) 



who suggested that the inner atmosphere of the sun is heated by acoustic waves that ar e 



generated in the convec tive zone. Later, theoretical and numerical studies (IStein 



1967 



Carlsson fc Stein 



1992) have demonstrated that acoustic waves are, in fact, abundantly 



generated in the convective zone and that these waves can, in principle, be responsible 
for chromospheric heating. However, measurements of the acoustic flux at different 
chromospheric levels usua 



ly fail to find sufficient energy to heat the whole chromosphere 
20051 ). However, as it has been shown recently, the chromosphere in 
the magnetic network may be heated by magnetoacoustic waves generated locally, inside or 



( iFossum fc Carlsson 



in the vicinity of the magnetic flux tubes ( IHasan fc van Ballegooijen 



20081 ). Also random 



Alfven waves can heat upper chromosphere via ion- neutral collisions and g enerate slow 



shocks, which can explain the formation of spicules (lErdelyi fc James 



200J). 



As an alternative explan ation for the chromospheric heating, it has been suggested 



( iParker 



1988 



Sturrock 



19991 ) that impulsive nano-flares, powered by magnetic reconnection 



events, could be responsible for chromospheric heating. Although the observations show 
numerous transient brightenings on the sun, these are insufficiently freq uent and insufficient 



y 



energ etic to explain the persistent UV emission of the chromosphere (lAschwanden et al 
20001 ). During solar flares, the chromosphere can be strongly heated and ionized locally 
by precipitating electron beams and evaporate upward, producing obse rved polarised Ha 



emiss ion via collisional interaction with neutral surrounding hydrogen ([Fletcher fc Brown 
1998k 
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Yet another possibility for 



currents ((Rabin &: Moore 



c hromospheric hea ting is the resistive dissipation of electric 



1984 



Goodman 



20041 ). Recent analysis of three dimensional 



vector currents and temperatures, deduced from spectropolarimetric observations of a 
sunspot from photospheric to chromospheric levels, has shown that, while resistive current 
dissipation can contri bute to heat the sunspot chromosphere, it is not the domina nt factor 



(So cas-Navarro 



1963 



Buneman 



2007). Recently, it has been suggested that the Farley-Buneman ((Farley. 



19631 ) inst ability (FBI), driven by convective motions, can be responsib le 



for chromospheric heating ((Liperovsky et al. 



2000 : 



Fontenla 



2005 



Fontenla et al. 



20081 ) 



The FBI is known to create plasma irregularities in the terrestrial ionospheric E-region, 
at heights where the electrons are strongly magnetized. The interplay of the earth's 
electric and geomagnetic field produces currents which give rise to the FB instability. 
Similarly, in those places where the electrons are strongly magnetized, the collisional drag 
of the ions by neutral flows can cause the dev elopment of a similar instability. Using the 



decrement of t 



he FB 



Fontenla et al. 



derived by ((Farley. 



19631 ) and assuming a negligible ion magnetization, 



(120081 ) concluded that the FB 



of the chromosphere. Earlier, the analysis of 



should be present at lea st in the upper half 



Liperovsky et al. 



(120001 ) had indicated that 



the FBI might operate in the chromosphere at heights h > 1000 km. 

However, the studies of the FBI in the chromosphere conditions are incomplete and 
they do not take into account two effects which under chromospheric conditions are 
important as we will show below. Firstly, if the finite magnetization of ions is taken into 
account, the Hall current perturbations weaken the FBI, and the system beco mes stable fo r 



any n eutral flow velocity when the ion magnetization factor k exceeds unity ( jFejer et al. 
19841 ) . Therefore, this instability can not operate in the upper solar chromosphere, where 
k > 1. Secondly, c ontrary to the E-layer plasma in the Earth's atmosphere, and top of the 



20071 ). the ionization degree in the solar chromosphere 



solar photosphere ((Petrovic et al. 
is quite high (10~ 2 — 10~ 4 ) and, consequently, Coulomb collisions can not be ignored as is 
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usually done in the study of the E-layer plasma. 

In this letter, we study the FBI taking into account both the finite magnetization of 
ions and the Coulomb collisions. We shall show that, in contrast to the situation in very 
weakly ionized plasmas, the relatively high degree of ionisation in the solar chromosphere 
makes the Coulomb collisions important for the FBI development. As a result, the 
instability becomes possible even in plasmas with an ion magnetization k larger than unity. 
However, by applying our analytical results to the solar chromosphere, we show that even 
though the FBI can sporadically appear in the chromosphere, it cannot be the main source 
of chromospheric heating. 



2. Formalism 

We consider a weakly ionized plasma consisting of electrons, one species of singly 
charged ions and neutral hydrogen. In the upper solar chromosphere, the positively charged 
particles are mainly protons, whereas at lower altitudes the positive charge is dominated 
by heavy ions. We therefore do not further specify the type of ions, so that our results are 
applicable to both the upper and the lower chromosphere. 

The dynamics of electrons and ions in such plasmas is governed by the continuity and 
the Euler equations, viz. 



dt 

and 



dn r + v • KV Q ) = o, (l) 



m a — = q a I E H 1 m e v ep {\ a - V Q * ) - m a u an (V a - V n ). (2) 

Here, a — e, i denotes electrons or ions, a* denotes the charged species opposite to a, at 
stands for e for the electron equation and for proton (p) for the ion equation. Also, n 
corresponds to neutrals, and n a denotes the density, V a is the averaged drift velocity, m a 
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is the mass, T a is the temperature, q a is the charge, u a p is the collision frequency, c is the 
speed of light, /C is The Boltzmann constant, and d a /dt denotes the convective derivative. 



For the electron-io n and electron-neutron collision frequencies we use the following 
1965) 



expressions (IBraginskii 



4(27r) 1 / 2 e 4 n e A 



3me /2 (/CT) 3 / 2 



(3) 



and 



O ~pr> 



<KT„ 



m r 



(4) 



where A corresponds to the Coulomb logarithm and in the former equation we take into 
account the fact that, regardless the mass of the dominant ion species, the ion-neutral 
collision frequency in the solar chromosphere is proportional to the thermal velocity of 
the neutral (hydrogen) component, because in the case of equal temperatures the thermal 
velocity of the neutral particles equals to or is greater than the thermal velocity of the ions, 
depending on the mass of the ions. The electron-neutral and ion-neutral collision cross 
sections are a Pn = 3.0 x 10 -15 cm 2 ( Bedersen fc Kieffer 



( jKrstic fc Schultz 



19711 ) and a in = 2.8 x 1(T 14 cm 2 



19991 ) . respectively. 



We assume, that the system is penetrated by a uniform magnetic field B and that the 
neutrals have a background velocity V n _L B. The Eqs. (j2J) then yield a stationary solution 
for the background ion drift velocity 

^1 + KKilpN) V n + K~V n X Z 



V,; 



(5) 



1 + 4 

and for the so-called 'current velocity', i.e. the relative velocity between the ions and the 
electrons, 

U = V i -V e = — — 5 . (6) 



Here k = uJcp/vpn is proton magnetization, 



OJr. 



1 + k( 

eB/m a c is the cyclotron frequency, 
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Ki = k(1 + ipN), ip = v en Vi n / uj cp uj ce and N = v ev jv en is the ratio of Coulomb and 
electron-neutral collision frequencies. 

On this background, we study linear electrostatic perturbations in the plane 
perpendicular to the background magnetic field. In order to simplify analysis we 
admit standard assumpt i ons for the study of t 



(jOppenheim et al. 



1996 



Schunk fe Nagy 



re FBI in the Earth's E-layer plasma 
2000). We assume quasi- neutrality (n e ~ n £ 



Technically, this means that, instead of using Poisson's equation, we use VJ = 0, where J 
denotes the electric current density. In addition, because n e and rii are indistinguishable, 
we use only one continuity equation. Finally, we treat the electrons as massless because the 
FBI occurs on an ion-neutral collision timescale which, for typical chromospheric conditions, 
strongly exceeds both the electron cyclotron gyration and the electron plasma oscillation 
timescales. 

We Linearize Eqs. ([I])-® and perform a Fourier transform of the obtained equations, 
and after long but straightforward algebra arrive at the following dispersion equation 



0. 



(7) 



where c s = [/C(T e + Tij/m^ 1 ' 2 is the sound velocity and u' in = m p Vi n /mi. In the derivation 
of the dispersion equation, we neglected all terms of the order of the small parameter 



if)K 2 ~ m e v en /m p i>p n ~ 2.6 x 10 3 . 

In the limit of low-frequency and long-wavelength perturbations (\u\ , |k • U | v' in ), 
we obtain the oscillation frequency uj t and the growth rate 7 of the Farley-Buneman type 
instability: 

= zr—y, (8) 
1 + ip 

and 



7 



K 2 /(l+N) 

(l + *p) 2 



k 2 c 2 s 
(k-C/ ) 



(9) 
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where ip = + N). 



3. Discussion and conclusions 

The terms proportional to iV = v ev jv en in Eqs. (EH) describe the effect of the 
electron-ion (Coulomb) collisions on the FBI and represent the main analytical result of 
this letter. When the Coul omb collisions are n eglected (N = 0), then the Eqs. (EH) reduce 



to the well known result of 



Fejer et al. 



(119841 ). which indicates that in plasmas with k > 1, 
the FBI cannot develop regardless of the neutral drag velocity. In contrast, if the Coulomb 
collisions are sufficiently frequent, the FBI can appear even when the ions are relatively 



highly magnetized. The depende nce of N on height, 



chromospheric model SRPM 306 (IFontenla et al. 



Dased on data of a semi-empirical 



20071 ). is shown in Fig. HJ It is seen that 



the Coulomb collisions dominate the electron- neutral collisions (N > 1) in the upper half of 
the solar chromosphere, at heights h > 1000 km. At these heights they cannot be ignored 
and facilitate the FBI. 

In the absence of Coulomb collisions, 



( lOppenheim et al. 



1996 



Schunk fc Nagy 



;he ph ysical background of the FBI is as follows 
200ol ): the electric field perturbation, 5E || k 
is the leading force acting on the ions and causes both the Hall (5Vf a11 ~ <5E x Bo) and 
the Pedersen (<5Vf ed ~ <5E) ion drift velocities. The FBI is powered by the force due to 
convective term (V j • V)(5Vf ed + <5V^ a11 ) in the ion momentum equation ([5]). This force, 
acting on the ions, causes the Hall and Pedersen responses in the perturbed ion velocity. 
The total velocity response parallel to k contains two parts: the Pedersen response (due 
to the Pedersen velocity <5Vf ed ) giving rise to the destabilizing term in ([9]), and the Hall 
response (due to the Hall velocity ^Vj 1311 ) giving rise to the stabilizing term proportional to 



K 
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Our analysis shows that, in the case of relatively frequent Coulomb collisions, the 
physics of the FBI is significantly modified. The reason for this is that the Hall and the 
Pedersen responses are both influenced by the Coulomb collisions but in a different way: 
the Coulomb collisions reduce the Pedersen response but not so much as the Hall response. 
This circumstance favors the FBI and makes it possible even for k > 1. 

We can determine the threshold value of the relative velocity necessary to trigger 
the FBI in the framework of the model SRPM 306. Fig. [2] shows the dependence of Uq t on 
height with (solid lines) and without (dashed lines) Coulomb collisions, for B = 30 G (thin 
lines) and for B = 60 G (thick lines). The left panel corresponds to the protons and the 
right panel to the ions with rrii = 30m p . 

The results obtained here allow us to draw conclusions about the possible role of the 
FBI in the inter-network chromospheric heating. The threshold value of the current velocity 
necessary to trigger the FBI corresponds to the current density Jo = eno e t/o r . Even for 
the lower chromosphere, where the positively charged particles are mainly heavy ions and 
Uq t ~ 2km/sec, very strong current densities Jq ~ 2.4 x 10 6 statampere/cm 2 are re quired 



Socas-Navarro 



(120071 ). at 



for the FBI to develop. According to recent observations of 
length-scales of the order 100 km and higher, the typical values of the observed currents 
are much smaller, ~ 5 x 10 4 statampere/cm 2 . It is in principle possible that such strong 
currents could exist at smaller scales. However, as we show below, in this case the heating 
rate produced by the frictional dissipation of the relative ion-neutral motion would be much 
higher than the power required to sustain the radiative loses in the chromospher e. Indeed, 



19651 ) 



the rate of frictional dissipation in partially ionized plasmas is ( iBraginskii 

Q fr = m e n e v ei (V e - V;) 2 + m e n e u en (V e - V^) 2 + m p n e v ei (V e - V;) 2 , (10) 
where the terms on the right-hand side are due to electron-ion, electron-neutral and 
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ion-neutral frictions, respectively. By substituting here Eqs. (JS])-([6]), we obtain 

Q fr = m e n e v en Ul (l + K 2 + N + ^ . (11) 

At relatively low heights in the solar chromosphere, where the positive charge is dominated 
by heavy ions, the friction dissipation is dominated by ion-neutral friction (the term 
proportional to 1/if))- Even for heights around h = 850 km, where Uq t ~ 2 km/sec 
and the electron density is relatively small in accordance to the SRPM 306, we have 
Qf r ~ 40 erg/cm 3 sec. In terms of the associated energy flux, such a heating rate would 
be produced by a wave flux F ~ Qf r H ~ 4 x 10 s erg/cm 2 sec dissipated in the lower 
chromosphere with H ~ 100 km being the characteristic width of the unstable layer. This 
value of the flux is at least one order of magnitude higher than necessary to compensate the 
radiative loses of the chromosphere and such fluxes are not observed in the chromosphere. It 
is not impossible, that the strong over-threshold currents occur sporadically at small scales, 
and drive sporadic FBI events at the length scales unresolvable for modern observations. 
However, given the strong frictional heating associated with such currents, the role of the 
FBI in the current dissipation and the associated heating is of minor importance. On the 
other hand, the small-scale plasma irregularities produced by the FBI can scatter radio 
waves, and hence provide a diagnostic tool for strong cross-field chromospheric currents if 
they exist. This last issue requires further investigation. 

In particular, our preliminary analysis of the full dispersion equation ([7]) has already 
demonstrated that for typical chromospheric parameters the strongest FB instability occurs 
at a characteristic wavelength that varies with height in the range 0.1-3 m. In the middle 
chromosphere at 1000 km height the maximum instability growth rate 7 ~ 2 x 10 3 s _1 occurs 
for waves with characteristic wavelengths A ~ 16 cm. The decimetric radio emission should 
effectively interact with the electron density irregularities produced by the FB instability, 
and this should result in observable scintillations of the decimetric radio emission. 
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In summary, we determined the threshold current velocity for the FBI to occur in a 
weakly ionized plasma taking into account the finite ion magnetization and the electron-ion 
Coulomb collisions. We have shown that, in the presence of Coulomb collisions, the FBI 
can occur even when the ion magnetization is greater than unity. Applying these analytical 
results to the solar chromosphere, we concluded that the FBI cannot be responsible for the 
chromospheric heating at global length scales. The FBI at small length scales cannot be 
excluded, but the heating produced by the FBI cannot compete with the frictional heating 
under chromospheric conditions. The small-scale irregularities generated by the FBI can be 
used for remote diagnostics of strong cross-field currents in the solar chromosphere. 
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Fig. 1. — The ratio of Coulomb and electron- neutral collision rates N as a function of height. 





400 600 800 1000 1200 1400 1600 
h, km 



400 600 800 1000 1200 1400 1600 
h, km 



Fig. 2. — Dependence of the threshold value of the velocity on height with (solid lines) 
and without (dashed lines) Coulomb collisions, for B = 30G (thin lines) and for B = 60G 
(thick lines). Left panel corresponds to the protons and right panel to ions with mi = 30m p . 



